Quantum systems in a stationary environment out of thermal equilibrium 



Bruno Bellomo, 1 ' 2 Riccardo Messina, 3 Didier Felbacq, 1 ' 2 and Mauro Antezza 1,2 

1 Universite Montpellier 2, Laboratoire Charles Coulomb UMR 5221 - F-34095, Montpellier, France 
2 CNRS, Laboratoire Charles Coulomb UMR 5221 - F-34095, Montpellier, France 
3 Laboratoire Charles Fabry, Institut d'Optique, CNRS, Universite Paris-Sud, 
Campus Poly technique, RD 128, 91127 Palaiseau cedex, France 
(Dated: January 16, 2013) 

We discuss how the thermalization of an elementary quantum system is modified when the system 
is placed in an environment out of thermal equilibrium. To this aim we provide a detailed investiga- 
tion of the dynamics of an atomic system placed close to a body of arbitrary geometry and dielectric 
permittivity, whose temperature Tm is different from that of the surrounding walls Tw- A suitable 
master equation for the general case of an iV-level atom is first derived and then specialized to the 
cases of a two- and three-level atom. Transition rates and steady states are explicitly expressed as 
a function of the scattering matrices of the body and become both qualitatively and quantitatively 
different from the case of radiation at thermal equilibrium. Out of equilibrium, the system steady 
state depends on the system-body distance, on the geometry of the body and on the interplay of all 
such parameters with the body optical resonances. While a two-level atom tends toward a thermal 
state, this is not the case already in the presence of three atomic levels. This peculiar behavior 
can be exploited, for example, to invert the populations ordering and to provide an efficient cooling 
mechanism for the internal state of the quantum system. We finally provide numerical studies and 
asymptotic expressions when the body is a slab of finite thickness. Our predictions can be relevant 
for a wide class of experimental configurations out of thermal equilibrium involving different physical 
realizations of two or three-level systems. 
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I. INTRODUCTION 

The absence of thermal equilibrium is a condition nat- 
urally verified in several biological systems as well as ex- 
perimental configurations [TH3]. Out of thermal equi- 
librium (OTE) systems have been recently subject to 
intensive investigations concerning heat transfer [3HT2] 
and Casimir-Lifshitz interaction pTH2Q] . It has been 
theoretically and experimentally demonstrated that sys- 
tems driven out of thermal equilibrium may show forces 
with new qualitative and quantitative behaviors: repul- 
sion, amplification of the force and stronger tunability 
[T3] [13] . Typical OTE configurations consist of an en- 
semble of bodies kept at fixed and different temperatures 
and embedded in a blackbody radiation emitted by the 
surrounding walls at a temperature different from those 
of the bodies, the whole system being in a stationary 
configuration. 

Recently, promising results have been also obtained re- 
garding the potentiality of OTE environments to control 
and manipulate the dynamics of atomic systems placed 
in this kind of environments [21]. This study may have 
experimental relevance for systems which are naturally 
out of thermal equilibrium, such as in recent studies 
concerning cold atoms close to superconducting surfaces 
[221 [23] , involving the tip of an atomic- force microscope 
AFM close to cold samples [24, 25] . in the case of atom- 
chip configurations [2HJ [21] and in all the setups aiming 
at investigating near- field properties [131 HU [28H3Q] . 

The lifetime 1/T of the atomic dynamics depends on 
the state of the local electromagnetic (EM) field with 



which the atom is coupled, i.e. on the fact that it is 
stationary or not, isotropic or not, that it is a vacuum 
state, a thermal state, a laser field, or other. Given an 
atom in a certain position, the local EM field is in general 
modified by the radiation emitted and/or scattered by 
surrounding bodies. This mechanism results in lifetimes 
depending on the relative distance between the atom and 
the bodies, on their geometrical and optical properties, as 
well as on their temperatures. Lifetime has been typically 
studied at thermal equilibrium, when the radiation emit- 
ted by the surrounding walls impinging the body is at 
thermal equilibrium with the body itself. Under this as- 
sumption several configurations have been investigated, 
among which are zero and non-zero temperatures, and 
infinitely thick planar and corrugated slabs [3TH33] . 

The atomic lifetime 1/T is strongly connected with the 
time evolution of atomic coherences, which may be nat- 
urally investigated by studying the reduced density ma- 
trix, also useful to evaluate the average values of phys- 
ical observables. Remarkably, it gives information on 
the atomic steady state and hence on the thermaliza- 
tion process. For instance, for configurations where the 
body is at thermal equilibrium with the environment at 
temperature T, the atomic density matrix evolves on a 
time scale 1/T to a diagonal matrix corresponding to 
a thermal steady state at temperature T. Contrarily to 
lifetime, and due to peculiar cancelations, at thermal 
equilibrium the atomic steady state depends only on the 
ratios hw nrn /k^T (where hw nrn are the energies of the 
internal transitions), being independent of the presence 
of the body. The richness of this system can be ex- 
alted if the atom is embedded in a stationary configu- 
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ration out of thermal equilibrium, made by a body whose 
temperature is kept fixed and different from that of the 
walls surrounding the body-atom system. The electro- 
magnetic structure holding in OTE systems has a com- 
plex nature, poorly investigated because of the impos- 
sibility of exploiting directly the fluctuation-dissipation 
theorem. Recently, using multiple-scattering theory and 
an appropriate use of the fluctuation-dissipation theorem 
at each different temperature present in the system, the 
radiation field in complex OTE configurations has been 
characterized in terms of its correlators [TTJ [12]. The 
knowledge of the correlators of the local EM field is the 
main ingredient needed to compute the transition rates 
regulating the atomic dynamics, as e.g. in the case of 
Kossakowski-Lindblad master equations in the Marko- 
vian regime [21] [34] . 

In this paper we provide a systematic derivation for 
the internal dynamics of an atom placed in front of an 
arbitrary body embedded in a thermal radiation whose 
temperature is different from that of the body. The pa- 
per is organized as follows. In Sec. [IT] we describe the 
physical model under investigation and we derive a mas- 
ter equation for the general case of a TV-level atom. In 



Sec. [Hi] we derive closed-form expressions for the tran- 
sition rates holding out of thermal equilibrium in terms 
of the scattering matrices of the body, valid for arbitrary 
geometrical and material properties. In Sees. IV and 
[V]we specialize our analysis to the case of a two- and a 
three-level atom, discussing how the atomic dynamics oc- 
curring at thermal equilibrium is modified by the absence 



of equilibrium. In Sec. [VTJwe derive explicit expressions 
for the transition rates for the case in which the body is a 
slab of finite thickness. This configuration is numerically 
investigated in Sec. |VII| where the cases in which slab is 
made of a dielectrics (GaAs) or a metal (gold) are con- 
sidered and compared. The dynamics of both two- and 
three-level atoms are discussed, showing the occurrence 
of peculiar phenomena emerging out of thermal equilib- 
rium, already pointed out in the general analysis. In Sec. 
I Villi we draw our conclusions. 



actual material constituting the walls. The atom inter- 
acts with the electromagnetic field (playing the role of 
bath B) generated by the walls and the body. The total 
Hamiltonian has the form 



H = H A + H E 



Hr 



(1) 



where H A in the free Hamiltonian of the atom and 
Hb the free Hamiltonian of the bath. We describe the 
interaction between the atom and the field using the 
multipolar coupling in dipole approximation [35] Hi = 
— D • E(R), where D is the atomic electric-dipole opera- 
tor and E(R) is the electric field at the atomic position 
R in the Schrodinger picture. In the interaction picture 
Hi becomes 



ffj(t) = -D(t).E(R,t), 



(2) 



where the time-dependent electric-dipole operator and 
electric field are defined by the transformations 
D(t) = exp(^H A t)T>exp(-^H A t) and E(R, t) = 
exp(^i?Bt)£ , (R) exp( — ^Hst). We describe the electric 
field in the interaction picture using a decomposition in 
which a mode of the field is identified by the frequency 
u;, the transverse wave vector k = (k x ,k y ), the polariza- 
tion index p (taking the values p = 1, 2 corresponding to 
TE and TM polarizations respectively), and the direc- 
tion or propagation = ±1 (shorthand notation <p = ±) 
along the z axis. In this approach, the total wavevector 
takes the form = (k, (j)k z ), where the z component 
of the wavevector k z is a dependent variable given by 

k z = \ ^ — & 2 , with k = |k|. The explicit expression of 



the field is 



E(R, t) = 2 Re 



r+oc 
JO 



duo 
2tt 



e -^*E(R,a;) 



(3) 



where a single-frequency component reads 

42\ 



E(R,c) = W ^c*^(k, W )^(k, W ), (4) 



II. MODEL 

We consider a TV-level atom A placed close to a body 
of arbitrary geometry and dielectric permittivity and em- 
bedded in an environmental radiation generated by the 
walls surrounding the system (see Fig. [TJ. The body and 
the surrounding walls have in general different tempera- 
tures, Tm and Tw, which are kept fixed in time realizing a 
stationary configuration out of thermal equilibrium. By 
assumption, the walls are far from both the body and the 
atom, their shape is irregular and they are at local ther- 
mal equilibrium at temperature Tw- As a result of these 
hypotheses, the radiation associated to the wall has no 
evanescent contribution reaching the body and the atom, 
and is a universal isotropic blackbody radiation. This 
blackbody radiation is by definition independent of the 
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Figure 1: (color online). The atom is kept fixed close to an ar- 
bitrary body whose temperature Tm is different from that of the 
surrounding walls, Tw- Tm and Tw are kept fixed in time, realizing 
a stationary environment for the atom. 
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E£(\l,lj) being the field amplitude operator associated 
to the mode (cj, k,p, <j>). For the polarization vectors ap- 
pearing in Q we adopt the following standard definitions 



e^ E (k,u;) = z x k = -(-k y x + k x y), 

e^ M (k,u;) = -e^ E (k,cj) x = - (-fez + # 2 k), 



(5) 



where x, y and z are the unit vectors along the three 
axes and k = k/fc. 



A. Derivation of the master equation 

Following [34] , we now derive a master equation for the 
reduced density matrix p{t) = TrB[ptot(£)] 5 where Tr# 
denotes the trace over the degrees of freedom associated 
to the bath. In the following we denote the eigenvalues of 
Ha by e and the projection onto the eigenspace belonging 
to the eigenvalue e by 11(e). The starting point of our 
derivation is, in the interaction picture, the von Neumann 
equation for the total density matrix ptot(^) : 



d 



Ptot(t) 



:[#/(*), PtotW], 



where iJj(t) of Eq. Q is rewritten as [34] 



(6) 



where i G {x,y,z} and Ai(u) = E e '-e=u; n ( e ) An(e'). 
Ai(uj) are eigenoperators of Ha belonging to the fre- 
quencies ±cj, being [iif^, A$(u;)] = — uAi(u) and 
= wAJH, and it holds A = £) w ^(cj) and 
A^(— cj) = The sum over uj appearing in ^ is 

extended over all the energy eigenvalues e' and e of Ha 
such that e' — e = cj. 

Using standard approximations, namely the Born, 
Markovian and rotating-wave approximations (see dis- 
cussion in Appendix [A]) , one obtains (using also the con- 
dition (^(r,^) = o)T 



dt 



pit) 



[H LS ,p(t)]+V(p(t)), 



(8) 



where the Hermitian operator (Lamb-shift Hamiltonian) 



H LS = S^Xi^'M^H^M 



(9) 



U l,J 



purely produces a shift of the atomic energy levels (we 
note that [Ha,Hls] = 0) and the dissipator term 
V(p(t)) has the form 



UJ l, J 



(10) 



where in Eqs. Q and ( [To| ) the frequencies uj assume 
both positive and negative values. In the two previous 
equations, jij(oo) and Sij(uS) are defined by 

i r°° 

^•M = — / dse wa (E i (B.,t)E j (B.,t-s)), 

Jo _ „ (11) 



From the previous equation it follows that S^-(cj) = 
+ ^ijM and 
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i r 



ds e iuJS (Ei (R, t)Ej (R, £ — s)). (12) 



The quantities are expressed in terms of the reser- 

voir correlations functions and in the case ^ is a sta- 
tionary state of the bath (as it is in our physical con- 
figuration), that is [Hb,Pb] — 0> they are homoge- 
nous in time. This means that (E^(R, t)i^(R, t — s)) = 
(E^(R, s)Ej(IL, 0)) and do not depend on time. 

Field correlators have been computed out of thermal 
equilibrium in a more general configuration, in pres- 
ence of a further object at a third different temperature 
[TTJ [12]. These depend on the three temperatures in- 
volved as well as on the material and geometrical prop- 
erties of the two bodies, taken into account by means of 
their classical scattering operators. We will show that by 
exploiting these results, the expressions of the transition 
rates can be explicitly given as a function of the reflection 
and transmission operators 1Z and T of the body M. 



B. A/-level atom 

Now we explicitly apply the previous derivation to 
an atomic system. We first consider the general case 
of an atom having TV states (ordered by increasing en- 
ergy) |1),|2),... and \N) : with corresponding frequen- 
cies cji, co>2, . • • , ojn- All the transitions between the N 
states are in principle allowed and the frequency differ- 
ence between two arbitrary levels n and m is indicated 
by uo nrn = uo n — uo m . We will assume from now on that n 
is always taken larger than m so that uo n m always repre- 
sents a positive frequency. The free Hamiltonian of the 
A^-level atom is 



N 



H A = ^2hu n \n)(n\, 



(13) 



n=l 



and the atomic dipole operator D(t) takes the form 



B(t) = ^(d mn \m)(n\e-^ t + d* mn \n)(m\e 

(14) 

where d mn = (ra|D|n) is the transition matrix element 
of the dipole operator between \m) and \n). In accor- 
dance with the previous remark, the sum runs over val- 
ues m, n G {1,2,..., N} such that n > m. Under these 



4 



assumptions, the sum over uj in Eqs. ([7]), Q and (10) 
runs over the values u nm and — uo nrn . By comparing Eqs. 
(§, Q and g§ one sees that (we remark that from now 
on > 0) 



where we have defined 



AM 



E 

{m,n}:a; 7^,7^=0; 



d mn |ra)(n| = A f (-cj), (15) 



meaning that the sum is over n and m such that u; nm = uj. 
In general several transitions can be associated to the 
same atomic frequency lj, both because of degeneracy 
and/or the occurrence of equidistant levels. By inserting 
Eq. ( [l5| in Eq. ([8| and coming back to the Schrodinger 
representation, one obtains the master equation in its 
most general form 



E w »h + E E E 

. n uj {m,n},{m f ,n'} i,j 



(sij(-uj) [d m / n /]i[d mn ]*<W \m) (m\ 

M[d |n)<n'|),p(t) 



7u(-w)[dm'n']i[d m „]* 



+E E E 

x (pmm> \n){n'\ - ^S nn t{\rnf)(m\,p(t)} 
+ 7ijM[d mn\i [dfTj/TT/jj 

x ^ n / n |m / )(m| - ^ mm /{|n)(n / |,p(t)} 



(16) 



where the sum X^{ m n} {m 7 n'} * s mean ^ over a ^ couples 
(m,n) and (m',n f ) such that cj nm = uj n > m > = w (be- 
ing alway n > m and ??/ > m'). The previous equation 
contains both terms coupling each transition with itself 
((m,n) = (m',n f )) and terms coupling different transi- 
tions having the same frequency ((m, n) ^ (m',n')). In 
the case of non-equidistant and non-degenerate levels, 
we have that to each transition corresponds a different 
value uj nm and we can write A(cj nm ) = d mn |ra)(n| = 
A^(— uJnm)- In this case, the previous master equation 
reduces to 



d 
dt 



pit) = -i ^y]cj n \n) (n\ +^S(-u nm )\m)(m\ 

n mm 

+ Y^S{u nm )\n)(n\,p{t)} 

m,n 

+ ^2r(-uj nm )(p mm \n)(n\ - ^{\m)(m\,p(t)}^J 
+ ^2 T (u nm )(p nn \m)(m\ - ^{|n)(ra|, /)(*)}), 



(17) 



0{ — ^nm) 


\ ^ / UJ 1 [J ] * 

— / j &ij y uJnm) [ u mn J i [ u mn J j 


17 I tlJn^ ) 




r(-w nTO ) 


— ^ ^ Tij ( ^nm) [dmn]i [dTTin] j 


r(^nm) 


= ^ ^ Tij (^nm) [dmn]j [dmnjj- 



(18) 



We now turn to the calculation of the functions 7^ (w) 
appearing in Eq. (18) and defined in (12). They are, 
using homogeneity in time, 

. . 1 f°° , f + °° du' f + °° dw" 

e- i ^"- u >{E i (R,Lj")Ej(K,cj')) 
+ e i("+»")°{El(R,Lo")E j (R,Lo')) 

~ ¥ J \(El(R, -uj)Ej(K,uj')) uj<0 

(19) 

where we have used (i^(R, uj)Ej(H, uj')) = 
(£|(R,cj)£j(R,u/)) = and ds exp(-ies) = 

2ir5(e). By using the decomposition of Eq. Q, we 
obtain for positive arguments of jij(oo) 



1 r°° du' ^ [ d 2 k [ d 2 \ 

xe^- K ' A )- R [^(k, W )],[<(k>')]* 
x(E$(k,oj)E^(k',oj')), 

whereas for negative arguments we have 



) 2 



(20) 



7*i(-^) 



x e 
x 



1 



+°° dw' 



E 



d 2 k f d 2 k' 



2tt ^ 7 (2tt) 2 7 (2tt) 2 
W ,p,p' J y ' J y ' 



^-K-).R^ (k;W)]|[ ^ (k>% 

« (k)W) ^; (k ' )W0)) 



(21) 



where we always assume uj > 0. We observe that last 
equation can be obtained by taking the complex conju- 
gate of Eq. ( 20 ) after having interchanged the operators 



E*(k,u) and^ t (k / ,o; / )- 



III. ATOM CLOSE TO AN ARBITRARY BODY 

At this point of the calculation, it is necessary to link 
the field correlation functions (£^(k, lj)E%, ^(k', uj')) and 
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(Ep '(k, cj)£y (k',u/)) to the two temperatures Tm and 
Tw and to the material and geometrical properties of 
the body as well. This can be done by following the 
derivation discussed in [12] in the more general case of 
two bodies and three temperatures. We are interested in 
the correlation functions of the field defined on the side of 
the body on which the atom is located (the right side, in 
our case). For a given set (cj,k,p), the expression of the 
the modes of the field propagating in the two directions 
is straightforward. We have, making the dependence on 
lj, k and p implicit, 



E~ = E^- 



(22) 



The field E~ propagating toward the body (i.e. toward 
the left) is only the environmental one, while the field 
E + propagating toward the right (see Fig. 0) results 
from the field T;( M )+ directly produced by the body, the 
transmission through M of the environmental field com- 
ing from the left, and the reflection by M of the one com- 
ing from the right. The operators 1Z and T are the stan- 
dard reflection and transmission scattering operators, ex- 
plicitly defined for example in [12] , associated in this case 
to the right side of the body. They connect any outgoing 
(reflected or transmitted) mode of the field to the entire 
set of incoming modes. 

Total field correlators are obtained by the knowledge 
of correlators of the fields emitted by each source to- 
gether with Eq. p2| ). The source fields have been char- 
acterized in [12] by treating each source independently 
as if it was at thermal equilibrium at its own temper- 
ature and thus applying the fluctuation-dissipation the- 
orem. The following symmetrized correlation functions 
((AB) sym = ((AB) + (BA))/2) have been obtained 



(4 M )+(k,c)4 M)+t (k^,^)) 



N(uj,T m ) 

ew) 



/sym " 2e c 2 
2ir5(cv - u')(p, k| (p^ ~ KP^P) + KP { °1 

P^tf -TP^tA |p',k'), 



(E^(k,u)E^'\k',u:')) sym = ^iV( w ,T w ) 



x 2nS(u - w')8m (p, k\P^ \p', k') 



(23) 



where 

N(.,T) = T coth(— 



+ n(to,T) , (24) 



with 



i(u,T) = (e^-l) \ 



(25) 



with II^ pw ^ and II ( ew ) being the projectors on the prop- 
agative (ck < cj, corresponding to a real k z ) and evanes- 
cent (ck > cj, corresponding to a purely imaginary k z ) 
sectors respectively. By combining Eqs. (22) and (23) 



we finally obtain the symmetrized correlation functions 
of the total field in the region of interest 



2e c 2 



(^+(k,^+ t (k / ,^))sym = 2^-^) 

(p, k| [n(u, T m ) (V^ ~ nV^rt + KV^ 
+ N(u>, T w ) (TV^T ] + KV^rt)] \p\ k'), 



2e n c 2 



(E+(k,Lj)E-i(k\uj')) sym = 2n6(LJ-Lj') 1 ^N(Lj,T w ) 

{E;(k,u:)E+\k',oj')) sym 

x (p,k|7>i7>^|p',k'), 
(E;(k,u)E-\k>,u')) sym 

x(p,k\P^\p>,k>). 

(27) 

To obtain the non-symmetrized versions of these cor- 
relation functions, appearing in Eqs. (20) and (21), 



2ir5(u>-u>')—^N(u>,T w ) 
2e c z 



2ir5(u;-u; f )—^N(u;,T w ) 
2e c z 



we first remark that the source correlation functions 
( [23] ) have been derived using thermal-equilibrium tech- 
niques at the temperature of each source individually 
(see [12 for a detailed discussion). As a consequence, 
we are allowed to use Kubo's prescription [36 , accord- 
ing to which in order to obtain (AB) from (AB) sym 
the replacement N(u,Ti) — > hu[l -f n(cj,Tj)] must be 
performed, whilst (BA) results from the replacement 
N(uj,Ti) — >• ftw n(uj,Ti). By exploiting the former re- 
placement rule, we obtain the modified version of Eq. 
d20b for uo > 



7ijM = 



uj 



rf—f 



d 2 k' 



s i(k-k')-r 



2e hc 2 ^ J (2tt) 2 J (2tt) 2 

^ k i| e ^-fc;*),[g+ (k)W)]i[ g+ (k>)] * 

"[1 + n(oj,T M )] (p^ ] ~ KV^tf 
[1 + n(u>, T w )] (tP^T* + KV^tf 



w) 



(pw) 



e+(k, W )],[e p -(k', W )]*^i P ; 



-i{k z -\-k*)z 



e;(k,o;)] i [eJ(k', W )]^T^ t 



-i{k z -k*)z\£ 



and 



r p {kMW-A^^)Y^-^]}\p'^)- 



(28) 

(p,k|P^ pw/ew) |p',k') = /c"(p,k|n (pw/ew) |y,k'), (26) We observe that jij(-cu) of Eq. ^ is obtained from 
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the previous equation by replacing [l + n(cj,T^)] with 
n(u,Ti) and by taking the complex conjugate. 

Using Eq. ( 28 ) one can finally cast the transition rates 
r(w nm ) and r(^co> nm ) of Eq. (18) under the form 



being in general T^ 171 ^ ^ T^ qj . No thermodynamical 
meaning different from the one above is associated to the 
notion of effective temperatures here introduced. The 
global dynamics can be then readily interpreted in terms 
of well-known thermal-equilibrium physics by means of 
the various ^ . However, out of thermal equilibrium, 
the various transitions feel in general different effective 
temperatures whose values depend on many factors such 
as the geometry of the body, the system-body distance 
and their interplay with the body optical resonances. By 
varying the various parameters one can tune the various 
effective temperatures far or close between them. 
In the case Tw 
n(uu. 



(pq) 



r(^nm) 

r(— cj nm ) 



where T (uj nrn ) 



nm J 



n 



(nm) N 

2 eff 
(nm) 

eff j 



(29) 



"L,|d ro „| 2 

37reohc 3 



is the vacuum spontaneous- 
emission rate relative to the transition \m) ^ \n). Here, 
we have introduced for an arbitrary transition frequency 



SlYC^^ [d mn ]i[dm n ]j f d 2 k f 

°wM = ^EE |dmn|2 

rt r>' i. n x ' 

i(k-k 



d 2 V 

(2tt)2 

(pw) 



P,P l :3 

; ')->,k|[ e -^-^^[e;(k, w )],[v(k', w )];^ 
<(*.+*:•)« [g + (k> w)] . [g - (k / ) w )]*7e P (pw) 

-e'( fe -^) z [6+(k, w )],[6+(k' )W )]* 



+ e 



a M 



2CJ ; — |u m)i | 

p,p' i,j 



I— I 

J (27T) 2 J 



d 2 V 



->(pw) 



: < ( k " k ') r <p, k| e**" )* [e+ (k, a;)], [e+ (k', a;)] * (p< 
- TiV^V) + PP^ - Pi^pt _ T p(pw) r t j | p ' ; k ') 



(30) 



and 



(nm) _ ^Km)^w)%Km)+wKm^M)«MKm) 



6 eff 



%Km) + «MKm) 



(31) 

The transition rates r(±co> nm ) of Eq. ( [29] ) are propor- 
tional to the term a^{uo nrn ) + ftMKm) which depends 
on material properties but not on temperatures, while 
^(nm) depends, in general, both on temperatures and the 
properties of the body. From the last equation it fol- 
lows that n(cj nm ,T min ) < n^ m) < n(cj nm , T max ), where 
^min = min(T M ,T w ) and T max max(T M ,T w ). As a 
consequence, the values of the transition rates r(±cj nm ) 
of Eq. p9] ) are always confined between their equilib- 
rium values at temperatures Tm and Tw- The transition 
rates relative to each frequency oj nm are equal to the 
ones at thermal equilibrium at an effective temperature 
T^ m \ associable to each transition. This is defined by 
n(cj nm ,T e ( ^ m) ) = n^ m) and is given by 



rp(nm) fc^nm 



ln(l + n<r )_1 ) 



(32) 



(nm) 



*eff 



T M = T, in Eqs. (|29j) and (|31|, 

,T) for any (nm). It follows that the 
transition rates r(±co> nm ) are factorized as a product of 
r (^nm), of a term (1 + n(uj nrn ,T) or n(uj nrn ,T)) de- 
pending on temperature, times the factor %(w nm ) + 
ftM(w nm ) depending on the properties of the body. In- 
deed, at thermal equilibrium the atomic decay rates de- 
pend on the presence of a material body because of the 
modification of the local field due to the field emitted by 
the body itself and to the way it scatters the one coming 
from the environment. 



Equations (17) and (29)-(31) allow one to investigate 



the time evolution of the atomic density matrix in the 
presence of an arbitrary body and for any couple of tem- 
peratures Tm and Tw- 

Remarkably, the explicit quantized form of Hb of Eq. 
has not been needed here. Indeed, all quantities we 
are interested in are related to the fluctuations of the 
fields, which under the local equilibrium assumption are 
provided, as explained in this section, by the fluctuation- 
dissipation theorem derived only by using general prop- 
erties of the fields. 

In order to get an insight on the atomic dynamics we 
consider in the next two sections two specific examples 
in which the atom has a simple level structure. 



IV. TWO-LEVEL SYSTEM 

In this section we assume that the atom has only two 
levels, ujo = UJ2 — oo\ being the transition frequency be- 
tween the excited state |e) = |2) and the ground state 
| fl > = |1) (see Fig. [fa)). 

The free Hamiltonian of the two-level atom is 



Ha = y^fiu n \n) M, 



(33) 



n=l 



and the atomic dipole operator D(t) assumes the simple 
form 



D(t) = d|l)<2|e- 



-iujot 



d*|2><l|e 



iujQt 



(34) 



where d = (g\T>\e) (we assume that its diagonal matrix 
elements vanish). In Eqs. ([9| and (10) the sum over i 
and j runs on z, j = x,y,z and the sum over uj on the 
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(a) | 2) 



|1) 



^32 




Figure 2: (color online). Scheme of levels and definition of 
transition frequencies for a two- (a) and three-level (b) atom. 



only two values cuo and —coo- By comparing Eqs. pj, ^ 
and (|34| one sees that A(uj ) = djl)(2| = A* (-Jo). By 
performing the sum over cj, Eq. (17) becomes 



j t p(t) = -i[£u Jn \n)(n\+S(-u Jo )\l)(l\ 

n 

+ S(wo)\2)(2\,p(tj\ +r(- W0 )(p 11 |2)<2|-i{|l)(l| ) p(t)}) 

+ r( Wo )(p22|i)<i|-^{|2)<2| )P (t)}), 

(35) 

where S f (— c^o), S(ojo), T(— cjo), r(co>o) are defined as 
in Eq. ( TTsb . The Lamb-shift Hamiltonian = 
5(— cjo)|1)\T[+ 5(cjo)|2)(2| induces a shift in the eigen- 
values of the free atom Hamiltonian Ha, which now be- 
come uji +S(—ujo) and 002 + S(ujq) with a difference given 
by A w = ujq + 5(cjo) — 5(— cjo)- r(cjo) and r(— ujq) are 
the transition rates associated to the down- and upward 
transitions respectively. From Eq. (35) the differential 
equations for p^ = (i\p(t)\j) follow: 



d 
dt 
d 

dt 
d_ 
dt 



Pll (t) = -r(-o;o)pii(*) + r(a; )P22(t), 



P22(t) = T(-(Jo)pil(t) ~ T(uj )p 2 2(t) 



(36) 



P12W. 



The solution of this equation reads 



Pn(t)\ _ /Wo) 



P22(t)J V^2(0) y 

Pi2(t) = pMoy^e 



1 



-7(w )* 



_ 7(wq) 



r(wd) 

7(w ) V r (~ w °) 

P2i(t) = p 12 (ty. 



(37) 



The decay rate 7(co>o) = r(co>o) + T(— cjo) is not influenced 
by the Lamb shift, which does not play any role also in 
the steady state of the system. At times t ^> 1/7(^0) the 
system thermalizes to the steady state 



p(t — > 00) 



1 



r(-^ ) + r(u;o) V r (-^: 



(38) 



The decoherence process is linked to the behavior of non- 
diagonal elements p\2 and p2i and is essentially regulated 



by the decay rate j{ljq)/2. As we will see, the behavior 
of the atomic evolution is qualitatively different between 
the equilibrium and non-equilibrium case. 



A. Case of thermal equilibrium 



When T M = T w = T, one has from Eq. (31) that 



(21) 

n eff = n (^o?^) so that the effective number of photons 
becomes independent of the material properties. More- 
over, the above stationary state (38) reduces to 



pit 



1 



l + 2n(w ,T) 



1 



n(to ,T) 
n(wo,T) 



(39) 

As a consequence of the factorization of thermal and ma- 
terial dependence in T(±uo), there is a simplification 
of the factor + ckm(^o) m the steady state. It 

follows that the asymptotic state is a thermal state of 
a two-level system at temperature T, independently of 
the presence and properties of body M. Out of thermal 
equilibrium, this peculiar cancellation does not hold any- 
more, the steady atomic state depending on the various 
parameters influencing the behavior of the T rates. We 
will show that the interplay between the body properties 
and the atomic position allows to realize a rich variety of 
asymptotic atomic configurations. 



B. Case out of thermal equilibrium 

When Tm ^ Tw, both the decay rates and the steady 
state of the system depend on the two temperatures and 
the properties of the body. In particular, these depen- 



dences are entangled in the term n^ ff of Eq. (31) (we 



(21) 



pose in the two-level case n eff 
in Eq. (38) we obtain 



nrff). Using Eq. (29) 



p(t 



1 



00 



1 + 2n e fi 



1 + n e ff 
n eff 



(40) 



We remark that the dependence on the material prop- 
erties has not disappeared, being indeed contained in 
the effective number of photons n e ff. Nevertheless, the 
steady state of Eq. (40) has the form of a thermal 



state so that the two-level system always thermalizes to- 
ward a thermal state even out of thermal equilibrium. 
The temperature T e ^ of this steady state is defined as 
in Eq. (32) by n(a;o,T e ff) = n e ff, it results equal to 
^B^eff = hoo/hi(l + l/n e ff) and is an intermediate tem- 
perature between Tm and Tw- 



We observe that the steady state (40 ) is defined by only 
one parameter (i.e. one of the two diagonal elements) 
being Tip = 1. As a consequence, any diagonal state 
satisfying the condition pn > P22 (always true in (40)) 
coincides indeed with a thermal state. In the case of a 
system characterized by more than two levels, a diagonal 
state is defined by more than one parameter and thus 
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can no longer be always identified with a thermal state, 
as we will see in the next section. 



are thus given by 



THREE-LEVEL SYSTEM 



We now focus our attention on the case of a three-level 
atom in a A configuration, previously considered in [21]. 
We label the three states with |1), |2) and |3), with fre- 
quencies co>i, co>2 and uo 3 (again in increasing order, see 
Fig. |2|b)). We assume that the interaction Hamilto- 
nian couples the states |1) and |3) as well as |2) and |3) 
(^32 7^ ^3i) 5 whereas the matrix element of D between 
|1) and |2) vanishes. The free Hamiltonian of the three- 
level atom is 



H A = y^fkj n \n)(n\, 



while the atomic dipole operator D(t) is given by 



D(f) = d 13 |l)(3|e-^ 3 ^ + d* 3 |3><l|e^ 3lt 
+ d 23 |2)(3|e— t + d5 3 |3)(2|e^ 2t , 



(41) 



(42) 



In Eqs. ([9]) and (10) the sum over uj runs now over the 
values 6J31, — 6J31, and — Co> 32 . By performing the sum 
over uo nm corresponding to the allowed transitions, Eq. 



(17) reduces to 



j f P= -*[5^Wn|n)(n 

n=l 


| + 5(-o;3i)|l)(l| 


+ (s(lo 31 ) + 5(W 3 2)) 


|3)(3|+5(-W32)|2)(2| )P 


+ r(-w 3 i)(pn|3)(3| 


- i{ll){ll,p}) 


+ r(w3i)(p33|l>(l|- 


|{|3><3|,p>) 


+ r(-w 32 )(p22|3)(3| 


-^{|2>(2|,p}) 


+ r(w 32 )(/933|2)(2|- 


|{|3> <3|,p>). 



(43) 



d , » 


-r( 


d . x 

Jt p22{t) = ' 


-r( 


|p3aW = r ( - 


-i 

d , . 

Jt^ = 


-r( 


iAj 


d , » 


zA £ 


^ 2sW = 


iA £ 



r(-cj 3 i) + r(-cj 32 ) 



r(cj 3 i) + r(-cj3i) + r(cj 32 ) 



*32 



r(a; 32 )+r(-cj 32 )+r(o; 3 i) 



P13W, 
P23(£). 

(44) 



The density matrix elements evolve asymptotically to 

r(u; 3 i)r(-cj 32 ) r(-cj 3 i)r(u; 32 ) 

pn(oo) = ~ , p 22 (oo) = ~ , 

Z Z 

r(-cj 3 i)r(-cj 32 ) . 

^33(00) = ~ , Pij(oo) =0 for z^j, 

Z 

z = t(uj 31 )t(-uj 32 ) + r(-cj 3 i)r(cj 32 ) 
+ r(-cj 3 i)r(-cj 32 ). 

(45) 

As far as the decoherence process is concerned, it is linked 
to the behavior of nondiagonal elements pi 2 , p\ 3 and p 23 , 
and then it is essentially regulated by the decay rates 
r(±cj 3 i) and T(±u 32 ). 



A. Case of thermal equilibrium 

In the case Tw = Tm 



(31) and noticing that = n(cj 3 i,T) and 



T, using Eqs. (|29) and 

f32) 



(73 32 ,T), the stationary state ( |45| ) reduces to 

pn (00) = n(o; 32 , T) [1 + n(u 31 , T)] /Z eq 

p 22 ((X)) = 77,(^31, T)[l + n(uj 3 2,T)]/Z eq 

^33(00) = n(cj 3 i,T)n(cj 32 ,T)/Z eq 

Z eq = 3n(cj 31 ,T)n(o; 32 ,T) +n(cj 31 ,T) + n(o; 32 ,T). 

(46) 

This state coincides indeed with the thermal state at tem- 
perature T of a three-level system, as expected indepen- 
dently of the properties of body M. 



The energy differences associated to the free Hamil- 
tonian modified by the Lamb-shift term read A31 = 
uj 31 +S(u 31 )-S(-uj 31 ) + S(uj 3 2), A 2 i = U2i + S(-u 3 2)- 
S(-uj 3 i) and A 32 = cj 32 + 5(o; 32 ) - S(-u 32 ) + S(u 3 i). 
The differential equations for the six diagonal and upper- 
diagonal matrix elements of the atomic density matrix 



B. Case out of thermal equilibrium 

The fundamental difference between the two-level and 
the three-level configuration is that in the latter we can 
associate an effective temperature comprised between 
Tw and Tm to each allowed transition, and that the 
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two effective temperatures do not in general coincide. 
Equivalently, to each allowed transition one can asso- 
ciate an effective number of photons given by Eq. (31) 



for (nm) G {(32), (31)}. The steady populations of Eq. 
( [45] ) can be rewritten in terms of the effective values of 
n associated to the two transitions as follows: 



0*11 (00) 

P22(oo) 
VP33(00), 

r = 3n 



1 

z 



,(31; 



''eff 



(1 



V 



(31)_(32) 
eff ™eff 



(31) 
n eff n eff 

,(31) 



(32) 



(47) 



<eff 



,(32) 
*eff ' 



The decay rates T(±uj nm ) have the form of Eq. ( |29| ) for 
(nm) G {(32), (31)}. Their values are always delimited 
by the equilibrium decay rates in the case Tw = Tm = 
T m i n and in the case Tw = Tm = T max . 



A first remark is that the steady state (47) is not in 



general a thermal state differently from the case of a 
two-level system. However a thermal state can be yet 
obtained if one succeeds in finding values of the various 
parameters such to make equal the two effective temper- 



atures, T 



(32) 



eff 



7 (32) 

reduces to (|46). 



^(31) 
L eff 



T. This gives n 



(31) __ 
eff — 



and n^ } = n(cj 32 ,T), so that the stationary state (47) 



A second relevant remark is that while the quantities 
T(±cj nm ), n^ 71 ^ and T^ m ^ associated to a given transi- 
tion (nm) are confined between their thermal-equilibrium 
values at T m i n and T max , this is not the case for the steady 
populations pn(oo) and p 22 (oo). It is possible to show 
that when n^ 2) = n(cj 32 ,T max ) and n^ 1} = n(cj 3i , T min ), 
Pi 1(00) reaches its maximum value which results to be 
larger than its value when Tw = Tm = T m i n (this is the 
maximum value for a thermal equilibrium configuration 
at a temperature comprised between T m i n and T max ). For 
this same condition p 22 (oo) reaches its minimum value 
which is smaller than its value when Tw = Tm = T m i n . 
This peculiar behavior is connected to the fact that the 
condition n^ 2) = n(cj 32 ,T max ) and nf^ = n(u 31 ,T min ) 
makes at the same time the transition |2) ^ |3) the 
most and |1) ^ |3) the less reactive. The opposite 

case is also of interest, namely when T e ^ 2 ^ = T m i n and 

T e ^ = T max . In this case the minimum of pu(oo) and 
the maximum of p 22 (oo) are obtained. Once again this 
result is respectively smaller and larger than their equi- 
librium values at Tw = Tm = T max . In this last case 



(31) 
if n(uJ 32 ,T u 

sion. 



n(u; 3 i,T max ) and n; ff = n(cj 32 , T min ), so that 
< n(co> 3 i, T max ) there is population inver- 



In the following we consider a specific example for 
which we give an explicit expression for «wKm) and 
aMK m ) providing then a numerical investigation of the 
peculiar effects described above for an arbitrary body. 



VI. ATOM IN FRONT OF A SLAB 



As a body we consider a slab of finite thickness 5 for 
which explicit expressions for the transmission and reflec- 
tion operators can be exploited [TTJ [12] . For this specific 
case, we will derive explicit analytic results for 
and chm(w). The atom has position R = (0,0, z) (we 
have chosen x = and y = in virtue of the cylindrical 
symmetry of the problem with respect to the axis z — 0) 
with z > 0, whereas the slab is defined by the two inter- 
faces z = and z = — S. As a result of the translational 
invariance of a planar slab with respect to the xy plane, 
its reflection and transmission operators, 1Z and T, are 
diagonal and given by 



(p,k|ft|j/,k'> = (27r) 2 5(k-k , )VPp(k^) ) 
(p,k|7>',k'> = (27r) 2 5(k-k')5 pp ,T p (k,uj). 



(48) 



They are defined in terms of the Fresnel reflection and 
transmission coefficients modified by the finite thickness 



r p (k,cj) 



2ik zm 5 

|(k, u)e 2ikzm6 



(49) 



r2(k,cj> 



,2ik z 



In these definitions we have introduced the z component 
of the K vector inside the medium, 



k2, 



(50) 



e(uj) being the dielectric permittivity of the slab, the or- 
dinary vacuum-medium Fresnel reflection coefficients 



k z -\- k zm 



e(uj)k z 



e(u)k z + k zm ' 



(51) 



as well as both the vacuum-medium (noted with t) and 
medium-vacuum (noted with t) transmission coefficients 



2k z 



k z H~ k zm 

2k zm 
k z -\- k zm 



2y/e(oj)k z 



e(uj)k z + 

2ygg^ 

e(uj)k z + k 



(52) 



After replacing the matrix elements (48) of the reflec- 
tion and transmission operators in Eq. (30) we obtain, 
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for an arbitrary transition frequency uo = uo n 



3c ^ d 




6+(k ) a;)] i [6+(k ) a;)] j (|p p (k,a;)| 2 + \r p (k,w)\ 2 ) 
[e+ (k, [e~ (k, w)],e 2ifc ^ / o J ,(k, u) 
[e p -(k, W )] i [£ p + (k, W )] J e- a ^p p (k,,)* 



«m(w) 



e p (k,w)]i[e p (k,w)] 

3c ^ [d m „]*[d 

mn 



Sttuj |d mn | 2 




x e 



x [e+(k,w)] i [e+(k,w)] i (l - |p P (k )W )| 2 - |r p (k, w )| 2 ) 

rlk k 

— [^(k )W )] i [e+(k,a,)]; 

2Im ^) z ( /0p (k, W )- /0p (k, W )*)|, 

(53) 

where we have used that fact that the polarization vectors 
are real quantities in the propagative sector (ck < ui). 
After performing explicitly the angular integrals, one ob- 
tains 



, , 3c ^ [d mn ]*[d mn ]j v-^ J [~< 



kdk 



x [M+(fc, 0/^(1 + |p p (fc, W )| 2 + |r p (fc )W )| 2 ) 
+ 2[M;(fc, W )] i Re( /0? ,(A : , W ) e 2ifc ^) 

3c ^ [d mn ]j[d m „]j „ x-^l H kdk 



Id I 2 



x [M^^l - Mfc, W )| 2 - |T p (fc, W )| 2 ) 

/ + °° A- rlk 
__ e -2 Im(*.)* [ M + (fc) w)] . Im(pp(fc) w)) 



where 



M?(fc,w) = (1,1,0), 



M£(fc,o;) 



|^| 2 ,0|fc z | 2 ,2fc 2 ). 



(54) 



(55) 



In a-w(oj nm ) and ftM(w nm ) only diagonal terms survive. 
They can be rewritten as 



, . l + B(cj nm ) + 2C(o; nm ) ~ 
«wKmj = " • a nm , 

/ v _ 1 - B(cj nm ) + 2D(cj nm ) ~ 

<^M(^nraJ — ~ • Cl nm , 



(56) 



where d„ 



= (\\d n 



1 = (1,1,1), we have used 



M[d nm y 2 ,|[d nm y 2 )/|d 



Ef c k dk , , /7 v 4cj 



2 

nm | 5 



(57) 



and introduced the vectors 



x (|p p (/c,cj)| 2 + |r p (/c,o;)| 2 ), 
3c ^ /~ /c dk. 



C ^ = ^z2 -r^M-(k,Lo)Re( Pp (k,to)e 



2ik z z\ 



P 

x Im(pp(fc,o;)) 



(58) 



The integrals appearing in B and C are defined in the 
propagative sector, while D is the only evanescent con- 
tribution. B is distance-independent and depends on the 
thickness 5, while C and D depend both on S and on the 
atom-slab distance z. Moreover, it is possible to show 
that C and D tend to for z —> oo, while for z —> 
C remains finite and D diverges (more precisely, its TM 
contribution) as (see Appendix [B] where the dependence 
of B(cj), C(uj) and D(cj) on S and z as well as the asymp- 
totic behavior of D(cj) are discussed) 



B(uj,z -> 0) 



3C3 Un(' U) ^] (1,1,2). (59) 



16uj 3 z 3 



s(uj) + 1 



It follows from Eq. (56) that for small z otM(u nm ) ^> 
%(w nm ) because of the divergence of D. As a conse- 
quence, all the effective temperatures tend to Tm in this 
limit (see Eqs. (31) and (32)) and thus the atom ther- 



malizes at the body temperature Tm- On the contrary, 
for large z, a^(uo nrn ) > «MKm) and the analysis of 



Eqs. (31 ) and (32 ) shows that the position of each 



(nm) 



in the interval [T^i n , T max ] is governed by B(cj nm ). From 
these comments it follows that it always exists a distance 
z for which %(w nm ) = ayii^nm)' At this position n^ 71 ^ 
does not depend on the geometry of the system and of the 
dielectric properties of the body (even if the point where 
this happens depends on them). This point delimits the 
two zones of influence where each temperature dominates 
for that specific transition. We also note that at thermal 
equilibrium B does not play any role, since all the quan- 
tities are proportional to the sum % + o;m, independent 
of B and asymptotically equal to 1 . We remark that since 
r(-w nm ) (see Eq. (29)) must be always larger than zero, 
by considering the cases Tw = Tm, Tw = 0, Tm ^ and 
Tw 7^ 0, Tm = three constraints follow (only two of 
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them are independent), for i = x : y, z: 



where z 



~ 2zuj 



. Equation (29) reduces to 



[B(w)]i - 

- [CHI* 



■2[C(4>0, 
2[D(w)]i > 0, 
f [D(w)]i > 0. 



(60) 



We note that the three relevant special cases in which 
the electric dipole is parallel (([d^]^ 2 = ([d^]^ 2 = 
|[d nm ]| 2 /2, [d nm ] z = 0) or perpendicular ([d nm ] x = 
[d nm }y 0, |[d nm ] 2 | 2 |[d nm ]| 2 ) to the surface, or 
isotropic (|[d nm ]i| 2 = |[d nm ]| 2 /3 for i = x,y,z), can be 
easily treated by specifying the symmetries of electric- 
dipole matrix elements in Eq. p0| ) or in Eq. ( [56] ) . 

Before providing a numerical investigation for both 
cases of two- and three-level atom, we briefly discuss 
some limiting cases already treated in literature. The 
case in which the body is absent is obtained by putting, 
in the general discussion of Sec. |III[ T = 1, 7Z = 0, or 
equivalently in Eq. ( [58] ) p p (k, uj) — and r p (k, uj) = 1 
(e(u) = 1), leading to B(J) = 1 and C(u) = D(cj) = 0. 
It follows that = 1 and qm(^) = 0, so that 

the rates r(±cj nm ) of Eq. (29) do not depend on 



the properties of the atomic dipole, and are equal to 
T((j) = T [1 + n(u;,T w )] and T(-uj) = T n(uj,T w ). As 
expected, the body temperature Tm disappears. 

In the limiting case of thermal equilibrium in the pres- 
ence of the slab, when Tw = Tm = T, Eq. (29) can be 
recast as 



r(o; nm ) 

r(— UJ nm ) 



roKm) 1 + C(o; nm ) + D(cj 

nm J 



l + n(w,T) 



n(w,T) 



(61) 



We notice that in the limit of zero temperature Eq. (61 ) 



gives the known-result for the modification of the atomic 
decay rate in front of a slab reported in [3TJ [32] . 

When the slab is a perfect mirror, one has e(u) —> oo, 
p p (k, uj) = ( — l) p and r p (k, uj) = 0. In this case the tem- 
perature of the slab does not play any role. This corre- 
sponds to the fact that under this assumption (dielectric 
permittivity real and infinite over the entire frequency 
spectrum) a temperature cannot be properly defined. By 
using in Eq. ( [58] ) p p (k,co>) = (— l) p and r p (k,cj) = 0, we 
have B(uj) = 1, ~D(uj) = and 



[C(w)]i 



3 r / sin z cos z \ , 
2[{— 



sin z , x 

[piX + iy ) 



5 iy + 2S i; 



nm / 

1 V ^nm ) J 

x ^1 + n(o; nm ,T w )\ 
n(cj nm ,T w ) y ' 



(63) 



Equations (62) and (63) allow one to retrieve known re- 
sults for this ideal case [31 . For example, for a dipole 
moment parallel to the surface the emission rate is re- 
duced to zero close to the slab, while in the case the 
dipole moment is perpendicular to the surface it is en- 
hanced with a factor 2 with respect to the absence of the 
slab. 



VII. NUMERICAL ANALYSIS 

In this section we provide numerical analysis for both 
cases of two- and three-level atoms treated respectively 
in Sees. [IV] and [V| considering an atom in front of a 
slab made of gallium arsenide (GaAs). The dielectric 
permittivity e{uS) of GaAs is described using a Drude- 
Lorentz model |38] 



e(uj) = e inf - 



■iTuj 



UJ* 



- UJ± 



iTuj' 



(64) 



(62) 



characterized by a resonance at u r = 0.506 x 10 14 rads -1 
and where ei n f = 11, uj\ = 0.550 x 10 14 rads _1 and 
T = 0.00452 x 10 14 rads _1 . This model implies a surface 
phonon-polariton resonance at uj p = 0.547 x 10 14 rads -1 . 
A relevant length scale in this case is c/uo r ~ 6 fim while 
a reference temperature is ftuj r /kB — 387 K. The main 
advantage of GaAs with respect to silicon carbide (SiC), 
previously used for the numerical analysis in [21] , is that 
its optically relevant frequencies are smaller than those 
of SiC. This allows one to observe interesting effects in 
a domain of frequencies where the thermal populations 
n(cj,T) are larger. Furthermore, in order to discuss the 
dependence of atomic dynamics with respect to changes 
of the overall optical response of the slab, we also briefly 
consider the case of a metallic slab. We specify here that 
all the numerical calculations presented in this section 
refer to the case of isotropic dipoles. 

As we will see, the atomic frequencies will be chosen 
of the same order of the slab resonances. In the case of a 
GaAs slab this can be achieved, for example, by consid- 
ering as an atomic system self-assembled quantum dots 
made of InGaAs, characterized by transition frequencies 
in the few terahertz range [39]. In these systems, the 
transition frequency can be modified by varying the size 
and the composition of the quantum dot. More in gen- 
eral, given an atomic system, one can look at materials 
for the body whose resonances match the atomic frequen- 
cies. 
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A. Transition rates and effective temperatures 

We first provide the analysis of quantities which are 
indeed common to the case of a two- and a three-level 
atom, namely the transition rates r(±cj) and the effective 
temperature for an arbitrary frequency uj. 

In Fig. [3] we plot T(-uj)/T (uj) and T(uj)/T (lj) as a 
function of the atom-slab distance z. For a given cou- 
ple of temperatures (T m i n , T max ) = (100, 600) K, we com- 
pare the two thermal-equilibrium configurations at 100 K 
and 600 K with the two possible configurations out of 
thermal-equilibrium. From the plot we see that out of 
thermal equilibrium the transition rates are always con- 
fined between the corresponding thermal-equilibrium val- 
ues. For small values of z only the temperature of the slab 
contributes. On the contrary, we see that the asymptotic 
value reached for large values of z differs from the equilib- 
rium realization at the environmental temperature. This 
means that even at large atom-slab separations both tem- 
peratures play a role in the atomic dynamics. 




10° 10 1 



Figure 3: (Color online) Main part: T(—uj)/Tq(uj) as a 
function of z for uj — 1.2cj r and 5 — 1 cm (semi-infinite 
slab) for (T m i n ,T max ) = (100, 600) K. Asymptotic curves for 
small z (black segments) for Tm — Tw — T max and Tm — 
Tw = T m i n . Inset: r(cj)/r (c<;) as a function of z. For 
both the main part and the inset: Tm = Tw = T max (red 
solid line), (Tw,Tm) = (T max ,T m i n ) (purple dot-dashed line), 
(Tw , Tm ) = (T m i 

n 5 T max ) (blue dotted line), Tm — Tw — T m i n 
(green dashed line). The symbols <(— indicate the asymptotic 
values (with respect to z) corresponding to the four couples 
of temperatures (see Appendix pi. 



Let us now discuss the dependence on the atomic tran- 
sition frequency uj for (T m i n , T max ) = (100, 600) K. Figure 
[4] shows that, both at and out of thermal equilibrium the 
rates T have a particularly pronounced dependence on uj 
around uj p . In this sense the atomic dynamics changes 
strongly if the transition frequency is tuned around the 
frequencies which are relevant for the dielectric response 
of the body (uj p in particular) . Figure [4] also shows that 
r(— uj) allows a much wider range of variations between 
the two thermal-equilibrium cases. This is due to the 
fact that, while Y(—uj) is directly proportional to n e fr, 




Figure 4: (color online). Panel (a): r(— uj)/To(uj r ) as a func- 
tion of uj. Panel (b): r(cj)/ro(c<;r) and Fo(uj) /Fo(uj r ) (gray 
small dotted line) as a function of uj. Same color convention 
of Fig. [3] for (T min ,T max ) = (100, 600) K. For both panels 
z = 1 /xm and S = 1.4 /mi. 



the dependence of T(uj) on 1 + n e ^ greatly reduces the 
relative changes with respect to Tm and Tw- 
in Fig. [5] we plot r(— uj)/To(uj r ) as a function of z 
and uj for S = 8.4 /xm and (T W ,T M ) = (600, 100) K. The 
plot elucidates the role of slab resonances uj r and uj p in 
the behavior of the transition rates. In particular, we 
observe that for any value of the atom-slab distance z, 
the transition rates always show a peak centered at uj p . 

In Fig. y the effective temperature T e ^ m ^ as a function 
of z and is plotted for four different frequencies. The 
temperatures chosen are again (Tw,^m) = (600, 100) K. 
The plot shows that by varying the atomic transition fre- 
quency uj nrn the effective temperatures may behave quite 
differently, with the common feature that for z small 
enough only the slab temperature remains relevant. For 
large z both temperatures contribute in a way which de- 
pends both on the frequency and the slab thickness 5. 
For intermediate distances the remaining range of tem- 
peratures is explored. We remark that small variations of 
frequency around the slab resonances (compare for exam- 
ple [6ja) and[6jb)) may lead to substantial changes in the 
behavior of 7^ m) . The comparison between the various 
plots may allow one to identify regions where the various 
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Figure 5: (color online). r(— uj)/Tq(uj) as a function of z and 
uj for 5 = 8.4 /im and (T W ,T M ) = (600, 100) K. 



transitions feel temperatures far between each other. 
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Figure 6: (color online). Density plot of T e ^ m ^ as a func- 
tion of z and S for four different frequencies for (Tw,Tm) = 
(600, 100) K. The chosen frequencies are: uo r (a), 1.02cj r (b), 
ujp (c), and 1.5 uo r (d). 



B. Two-level atom 

We now focus our attention on the case of a two-level 
atom. For this configuration, the steady state of the sys- 
tem is of course identified by only one parameter, for 
example the excited-state population ^22(00)- In Fig- 
we plot this population as a function of z and uoq for 
8 = 1cm (semi-infinite slab). In the main part, the cho- 
sen temperatures are (Tw,^m) = (200, 500) K. On the 
left (right) of the main part, the density plot at thermal 
equilibrium at 500 K (200 K) is given, which is indepen- 
dent of z. The main part shows that out of thermal equi- 
librium P22 depends on the atomic position. From the 
comparison of the three plots, we confirm that for small 
z only the slab temperature contributes, while for large 
z both Tw and Tm are relevant to determine the values 
of the population. The dashed lines (corresponding to 
uo r and Up) highlight that the variations of the atomic 
population are paricularly pronounced around the slab 
resonances. 
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Figure 7: (color online). Density plot of the excited-state 
population ^22(00) as a function of z and cjo for 5 = 1cm, 
corresponding to a semi- infinite slab. The temperatures are 
(T w , T M ) = (200, 500) K. On the left (right) of the main part, 
the density plot at thermal equilibrium at 500 K (200 K) is 
given. 

In Fig. |8]p22(oo) is plotted as a function of z for two 
values of 6 at ujq — 0.9 uj r for (T m i n , T max ) = (200, 500) K. 
We see that for S = 0.01 jam the atom thermalizes at large 
distance z at a temperature very close to Tm (see Fig. 
|8ja)), while by increasing the thickness to 300 /im (limit 
of semi-infinite slab) one obtains a steady state at an 
intermediate temperature between Tm and Tw (see Fig. 
|8jb)). This shows that the slab thickness S plays a role in 
defining the region of influence of the slab temperature. 

The mathematical reason for which in close proxim- 
ity of the body the atom tends to thermalize to the 
temperature of the body is the divergence of the inte- 
gral D(cj) of Eq. (58), corresponding to the fact that 



the field is dominated by the evanescent radiation pro- 
duced by the body and confined near its surface. The 
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Figure 8: (color online). £22(00) as a function of z for cjo = 
0.9o; r , (T m in,T max ) = (200, 500) K. Same color convention of 
Fig. [3] The two considered thicknesses are S = 0.01 /im (a) 
and S — 300 /im (b). 

absence of an evanescent contribution associated to the 
surrounding walls results in an atomic state which, even 
for large atom-body distances, does not tend to a thermal 
state at temperature Tw. Far from the surface, the atom 
sees only propagative contributions to the radiation (in 
other words, the atom is always far from the surrounding 
walls), and its state remains dependent on the tempera- 
ture and properties of the body. 

C. Three-level atom 

We now turn our attention to the case of a three-level 
system. As we will see, in this configuration the presence 
of two different effective temperatures (each one associ- 
ated to one of the two allowed atomic transitions) allows 
to produce a wider variety of physical situations, and 
in particular the population inversion of the two lowest- 
energy states. 

In Fig. [9ja)-(b) we plot the steady populations as 
a function of z for two different values of <5, choosing 
(Tyv, Tm) = (300, 50) K, u 32 = 1.02 u r and u 31 = u p . 
The plots show the occurrence of population inversion 
between the states |1) and |2). It is also shown that by 
varying only the slab thickness, for values of z around 
0.5 /im one can pass from a steady state with almost 
all the population in |2) to another one with almost all 
the population in |1). This kind of behavior corresponds 
to high values of purity which quantifies how pure the 
atomic state is [34] . This is denned for an arbitrary three- 
level system by 

U(p) = Tr[p 2 ] = /4 + 2(M 2 + \pis\ 2 + \P23\ 2 )- (65) 

i 

In the steady configuration the state is diagonal and high 
values of purity are reached when pn(oo) or /?22(oo) tends 
to 1 (assuming that ^33 (00) remains much smaller than 
both pn(oo) and £22(00)). In Fig. [9jc)-(d) the density 
plot of purity is represented by varying the two frequen- 
cies CJ32 and CJ31 for z = 0.47 /im and 5 = 0.01 jam in (c) 
and S = 2/im in (d). The black dotted lines correspond 
to uo r and uo pi pointing out that high values of purity can 
be obtained only when at least one frequency is close to 
slab resonances. 



z (jum) z (urn) 




1 1 iii- 1 1 : 1 

0.6 0.8 1 1.2 1.4 0.6 0.8 1 1.2 1.4 

o) 32 / oo r (ti 32 l co r 



Figure 9: (color online). Steady populations (red solid line 
for pn(oo), green dashed line for ^22(00), purple dot-dashed 
line for /?3 3 (oo)) as a function of z for 5 = 0.01 /im (a) and 
5 = 2/im (b), being (T W ,T M ) = (300, 50) K, cj 3 2 = 1.02 cj r 
and CJ31 = Up. Density plot of purity as a function of the two 
frequencies CJ32 and CJ31 for z = 0.47 /xm, S = 0.01 /im (c) and 
S — 2/im (d). The black dotted lines correspond to uo r and 
ujp and highlight the zones where high values of purity are 
obtainable. 



In order to stress the occurrence of inversion of pop- 
ulation ordering and its dependence on the atom-slab 
distance and the slab thickness, in Fig. [To] we present 
a density plot of the ratio P22(oo)/pn(oo) as a function 
of z and S for (T W ,T M ) = (500, 100) K, u; 3 2 = 1.02 cj r 
and co>3i = 1.57 uj r . The dotted line corresponds to a 
ratio equal to 1, that is when the two populations co- 
incide. This plot points out a strong dependence of the 
ratio on the two involved parameters, with regions where 
the ratio is larger than one characterized by inversion of 
population. This indeed occurs even if one of the two 
frequencies is far from uo r . 

In Fig. pT[a) we trace the surface of points 
(pn(oo), P22(coJ) by varying z in [0.01, 100] /im and S in 
[0.01/im,lcm] for (T W ,T M ) = (600, 50) K, U32 = 1.02 cj r 
and co>3i = ujp. The region covered by the steady states 
is compared with the curve of thermal states between 
Tw and Tm (red continuous line) and outside them (yel- 
low dashed curve). The plot proves that a large va- 
riety of steady states can be obtained. In particular, 
for S > 3 /im (semi-infinite slab) and z = 0.54 /im, 
the steady state is almost a thermal state at T ~ 5K 
(this value is obtained by looking for the closest ther- 
mal state by using the trace- norm distance). This cool- 
ing mechanism is connected to the behavior of the effec- 
tive temperatures, which for these values of z and S are 
T e ( ^ 2) ~ 145 K and T^ 1] ~ 59 K. We remark here that 
the cooling concerns only the internal atomic degrees of 
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Figure 10: (color online). Density plot of p22(oo)/pn(oo) 
as a function z and 6 for (Tw,Tm) = (500, 100) K, CJ32 = 
1.02 uj r and = 1.57cj r . The dotted line, corresponding 
to p22(oo)/pn(oo) = 1, highlights a large regions of occur- 
rence of inversion of population ordering as well as a strong 
dependence on z and S. 



freedom and not the external motion. In Fig. [TTJb) the 
same plot is presented by inverting the two temperatures 
(T W ,T M ) = (50, 600) K. The resulting shape of the sur- 
face is completely different, showing that the behavior 
of the steady state is not symmetric for exchange of the 
two temperatures. In particular, for large values of z the 
steady state remains quite far from the thermal state at 
T w . 
In Fig. 
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we show that, by varying only Tw and 
Tm in the interval [50, 600] K while keeping all the other 
parameters fixed, large regions in the space of states can 
be explored. Different shapes are realized when one varies 
for example usi and CJ32, as shown in the inset. Figure 
12 confirms that non-equilibrium configurations provide 
new tools to realize a large variety of non-trivial steady 
states. 



D. Gold slab 

In this subsection we briefly discuss the qualitative 
modifications to the atomic dynamics induced by the re- 
placement of GaAs with a metallic surface. To this aim 
we consider a slab made of gold, whose dielectric permit- 
tivity e(uj) is described using a Drude model 



e(w) = 1 



(66) 
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Figure 11: (color online). Region occupied by the couples 
(pn(oo), £22(00)), for z varying in [0.01, 100] /am and S in 
[0.01 /zm,l cm], for CJ32 = 1.02o; r and CJ31 = uj p . The red 
curve represents the thermal states between Tw and Tm while 
the yellow dashed curve corresponds to temperatures exceed- 
ing the ones involved in the system. For each value of 5 we 
represent a continuous curve by varying z. In (a) the chosen 
temperatures are (Tw, Tm) = (600, 50) K while in (b) we have 
(T W ,T M ) = (50, 600) K. 



model implies a surface plasmon resonance at uj p = 
96.987 x 10 14 rads -1 . Bothcj p i andco> p correspond to tem- 
peratures much higher than a typical room temperature 
of 300 K. However, in order to make effects out of ther- 
mal equilibrium emerge, we prefer to limit our analysis 
at frequencies around room temperature, that is around 
cj r = 300k B /h ~ 0.392 x 10 14 rads -1 . In this range of 
frequencies we observe much less richness in the behavior 
of the quantities of interest as a function of the various 
parameters. 

Analogously with what presented in Fig. |6j we plot 



characterized by a plasma frequency uj p \ = 137.2 x 
10 14 rads _1 and where T = 0.4059 x 10 14 rads -1 . This 



in Fig. 
of z and" 
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the effective temperature T^ 171 ^ as a function 
for four different frequencies and (Tw,2m) = 
(600, 100) K. By varying u nm around uj r the only visible 
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Figure 12: Points (pii(oo),p22(oo)) for all possible couples of tem- 
peratures (T W ,T M ) in the interval [50K,600K] for S = 300 /im 
(semi-inifinite slab) and z = 1 /im. Thermal states are repre- 
sented by red solid line (50 K < T < 600 K) and yellow dashed 
line (T < 50 K and T > 600 K). In the main plot, uu 3 2 = 1-02 uj r 
and = uj p . In the inset, CJ32 = 0.9 uj r and = uj r - Darker 
zones correspond to a higher density of states. By moving along 
the border of the surface we have the following variations of tem- 
perature. From the comparison with the inset one sees that in 
correspondence with the same change of temperatures, different 
trajectories are obtained. From 'a' ((Tw, Tm) = (600, 600) K) to 'c' 
((T w , T M ) = (600, 50) K), T w is fixed while T M goes from 600 K to 
50 K. From 'c' to 'b' ((T w , T M ) = (50, 50) K), T M is fixed while T w 
goes from 600 K to 50 K. From 'b' to 'd' ((T W ,T M ) = (50, 600) K), 
Tw is fixed while Tm goes from 50 K to 600 K. From 'd' to 'a', 
T M is fixed while T w goes from 50 K to 600 K. 'e' ((T W ,T M ) = 
(250, 200) K) represents a point internal to the surface. 



effect is a global horizontal shift of the plot, correspond- 
ing to the fact that by increasing the frequency the region 
where the effect of the slab dominates decreases. This ef- 
fect is connected to the rapid decay with lj of the imag- 
inary part of the dielectric permittivity, which in turn 
reduces the contribution of the term D(co> nm ) (the only 
one diverging at smaller distances) to ftM^nm) (see Eq. 
( [56] ) and Appendix [5) . 

In the case of a two-level system, in analogy with Fig. 
[7j we plot in Fig. 14 the excited-state population ^22(00) 
as a function of z and ujq for 5 = 1 cm (semi-infinite slab). 
The plot highlights a smoother dependence of ^22(00) on 
z. Moreover, we see that in the case of a metal thermal- 
equilibrium values at Tw are always recovered for large 
z. This is technically due to the fact that in Eq. (56) the 
contribution ~B(u) n 



i) has always values close to 1 so that 
in the zone where it dominates (that is asymptotically in 
z), it makes only %(w nm ) contribute and thus the slab 
temperature irrelevant. 

In the case of a three-level system, in analogy with 
Fig. 10 , we plot in Fig. [l5] the ratio p22(oo)/pn(oo) 
as a function z and S for (Tw,Tm) = (600, 60) K, CJ32 = 
0.704 cjr and CJ31 = co>r. Black dotted lines indicate where 
the ratio is equal to 1, that is when the two populations 
coincide. The comparison between this figure and Fig. 
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Figure 13: (color online). Atom in front of a gold slab. Den- 
sity plot of T^ m ^ as a function of z and 8 for four different 
frequencies for (Tw,Tm) = (600, 100) K. The chosen frequen- 
cies are: 0.1 cjr (a), 0.5 uj-r (b), cjr (c), and 2.5 cjr (d). 








0.05 


1 0.15 


2 0.25 


3 0.35 


4 0.45 



4 

3h 




10 u 

z (|im) 



Figure 14: (color online). Atom in front of a gold slab. Den- 
sity plot of the excited state population ^22(00) as a function 
of z and ujo for 5 — 1cm (semi-infinite slab). The tempera- 
tures are (T W ,T M ) = (200, 500) K. On the left (right) of the 
main part, the density plot at thermal equilibrium at 500 K 
(200 K) is given. 



10] shows that in this case the dependence of the ratio on 
the two involved parameters is less rich and inversion of 
population is reached in a small region of the (z, S) plane, 
with ^22(00) only slightly larger than pn(oo). 
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Figure 15: (color online). Atom in front of a gold slab. 
Density plot of ^22(00) / pn(oo) as a function z and 6 
for(T w ,T M ) (600, 60) K, cj 3 2 0.704 cj r and cj 3 i u; R . We 
observe a weak dependence on z and with small regions of 
occurrence of inversion of population ordering. 



have provided a detailed numerical investigation for both 
the cases when the slab is made of a dielectric and a 
metal, pointing out that dielectric configurations present 
a richer variety in controlling the atomic dynamics. The 
slab thickness regulates the distances from the slab at 
which the evanescent contribution of the field emitted by 
the slab dominates the atomic dynamics which, in this 
limit, becomes similar to an equilibrium dynamics at Tm. 
Moving the atom far enough from the slab, the effect of 
Tw becomes relevant. This effect depends on the relative 
value of the frequency of the transition with respect to 
the slab resonances. As a consequence, for each transition 
the influence of the two temperatures can be also quite 
different, allowing a strong manipulation of the atomic 
dynamics by varying Tm and Tw- 

The results reported in this paper could be of interest 
for experimental investigations in the absence of ther- 
mal equilibrium involving real or artificial atoms (such 
as quantum dots) . Our work suggests that similar effects 
will be present also in the case of more atoms trapped 
close to a substrate [14] or traveling in an atomic beam 
[40] in proximity of a surface. For quantum dots one can 
also consider mechanical manipulation of the dot [41]. 
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VIII. CONCLUSIONS 

In this paper we have reported a detailed investigation 
of the dynamics of an elementary quantum system placed 
close to an arbitrary body whose temperature Tm is kept 
fixed and different from that of the surrounding walls 
Tw- We have derived a suitable master equation for the 
atomic dynamics in such stationary environments out of 
thermal equilibrium, providing closed-form analytic ex- 
pressions for the transition rates governing the dynamics 
in terms of scattering matrices of the body. We have 
pointed out relevant differences with respect to the case 
of thermal equilibrium, with the steady state now depen- 
dent on a complex interplay between the atom-body dis- 
tance and the geometrical and material properties of the 
body. The overall dynamics can be readily interpreted in 
terms of effective temperatures associable to each allowed 
transition. After treating the general case of a TV-levels 
atom, the case of two- and there-level atoms have been 
discussed in detail. In the first case, the atom thermalizes 
to a thermal state at a temperature comprised between 
Tm and Tw, while in the second case the steady state 
is not thermal in general and the steady population may 
exceed the equilibrium values at Tm and Tw- This effect, 
emerging out of equilibrium, allows one to strongly ma- 
nipulate the steady state, realizing inversion of popula- 
tion ordering and cooling of internal atomic temperature 
at values external to the interval between the involved 
temperatures. We have then specialized our analysis to 
the case in which the body is a slab of finite thickness, 
deriving explicit expressions for the transition rates. We 
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Appendix A: Approximations 

In this appendix we briefly discu ss the approximations 

EH. 



used in the derivation of Sec. II A 



The first approximation consists of a weak-coupling 
assumption which allows to expand the exact equation 
of motion for the density matrix up to second order 
in the coupling constant. Together with the condition 
Ptot(t) — p(t)pB this leads to the Born approximation 
of the master equation. The second approximation is 
the Markov approximation in which the quantum master 
equation is made local in time: the density matrix p(t), 
in principle depending on p(s) for any s < t, is assumed 
to depend only on the density matrix at the same time t. 
Furthermore, the master equation is made independent 
of the initial state of the system. This Born-Markov limit 
is valid if the bath correlation time tb is small compared 
to the relaxation time of the system tr, that is tb «t^ 
Finally, in the rotating wave approximation, rapidly os- 
cillating terms proportional to exp[i(a/ — uj)t] for uo' ^ uo 
are neglected, ensuring that the quantum master equa- 
tion is in Lindblad form. The corresponding condition 
is that the inverse frequency differences involved in the 
problem are small compared to the relaxation time of 



the system, that is ta 



<C tr. The master 
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equation given in Eq. ([8| can be used only for values 
of the various parameters (coupling constants, tempera- 
tures, atom frequencies, etc.) such that the conditions 
stated above are satisfied. 

We first estimate tr for a given transition of frequency 
u nm . For a two- level atom, we argue from Eq. (37) that 
it is given by 

tr ~ [iKm)] 1 [r(-o; nm ) +T(u; nm )] \ (Al) 

For an atom with more than two levels, tr will be a 
complex function of all the various transition rates. For 
our numerical estimation, we consider the typical case 
of dipole momentum of the order of 10 -29 C m. In or- 
der to test the validity of the Born-Markov approxima- 
tion, we focus on the worst case. The smallest values 
of tr are obtained for uj nm = uo pi for small values of 
z (which in our numerical analysis start from 0.01 /im) 
and for large values of S. Let us fix the temperatures as 
T max = 600 K, T min = 100 K and S =1 cm. The smallest 
value of tr, obtained for z = 0.01 /im and Tm = T max , 
is tr ~ 4 x 10 -11 s. In this limit, the results at equilib- 
rium (Tw = T max ) and out of equilibrium (Tyy = T m [ n ) 
are equal (only the slab temperature is relevant). By 
increasing z, tr rapidly increases. At z ~ 0.04 pm we 
have tr ~ 2 timesl0~ 9 s. Also by slightly moving uj nrn 
around uo p we observe a rapid increase of tr. For example 
at uj n m — 1-04 uj p and z ~ 0.01 /im, r^~4x 10 -9 s. At 
Unm = Up and large values of z, we find very large times 
being tr ~ 0.005 s for Tw = T max . Similar considera- 
tions and numerical values are found for the three-level 
case in a A configuration analyzed along the paper. 

We now estimate tb by looking at the decay time 
of (Ei(R,,s)Ej(H,0)). This means computing (compare 
with Eq. (pb) 



(E l (R,s)E J (R,0)) 



2tt 



/+oo 
-oo 



s 7i ». (A2) 



We can thus numerically estimate t# and compare it with 
tr. In absence of the body the decay time tb is fixed 
by the thermal fluctuations giving times for tb of the 
order of 10 -13 s. In presence of the body, tb has a weak 
dependence on z and the two temperatures, and assumes 
values of the order of (1 — 10) x 10 _12 s. These values 
are almost always much smaller than r#, becoming close 
only for z ~ 0.01 /im. Once again, close to the slab only 
the slab temperature is relevant, so that the analyses at 
equilibrium and out of equilibrium give analogous results. 
Being the condition tb <C tr always satisfied for values 
of the dipole momentum small enough, then the Born- 
Markov approximation is well satisfied in our system. 

With regards to the rotating wave approximation, in 
the case of a two- level system ta ~ (2cj nm ) _1 , so that one 
has to check if (2cj nm ) _1 <C [T(-uj nrn ) + T(uj nm )] . In 
our system, the frequencies investigated are of the order 
of (0.05 — 1.5) x 10 14 rads _1 , so that ta is of the order 
of (0.3 — 10) x 10 _14 s. These values are typically much 
smaller than tr, which, as discussed before, in the worst 



case is of the order of 4 x 10 11 s. Similar considerations 
hold for the three-level atom considered in the paper. 



Appendix B: Integrals B(cj), C(u) and D( 



This section is devoted to a brief discussion of the de- 



pendence of B(co>), C(cj) and D(cj) (defined in Eq. (58)) 
on the atom-slab distance z and the slab thickness S. 

Let us start with the contribution B(cj), the only one 
appearing both in aw and As remarked before, it 
does not depend on z but only on the thickness S. Using 
Eq. ( |57| ) and noting that the sum \p p (k, uj)\ 2 + |r p (fc, uj)\ 2 
belongs to [0,1], we have as an immediate result that 
< [B(o;)]i < 1 for i = x,y,z. Moreover, for S — >• +oo 
we have T p (k,u) — > and we can replace p p (k,u) with 
the bulk Fresnel coefficient r p (k,u), obtaining 

B H = £E jf ^M+(fc, w )Mfc,c)| 2 . (Bl) 



Finally, for 5 = (namely in the absence of the slab) we 
have p p (k,u) = and T p (k,u) — 1 so that [B(cj)]i = 1. 
We observe that in the same limit we have G(uj) = and 
D(cj) = 0, so that = 0, i.e. the dependence on 

the properties of the slab naturally disappears. 

Let us now discuss the properties of C(cj). We first 
observe that C(u) is finite in the limit z — >• 0, as a conse- 
quence of Lebesgue dominated convergence theorem. As 
for the limit z — > +oo, we set k = — , z = — and we 

make the change of variable y/l — k 2 — » s. The integral 
C(lj) is reduced to the form 



](s)e 2is ~ z ds, 



(B2) 



for some smooth function g(s). An integration by parts 
then shows that C(u) — ^ as z —> +oo. 

The asymptotic behavior of D(cj) for z — >> +oo is ob- 
tained by noting that the integrals involved are, after the 
change of variable — 1 — » 5, of the form 



r + oo 

Jo ' 



ds, 



(B3) 



where the function g(s) has a polynomial behavior near 
5 = 0, i.e. g(s) ~ Q s n for some non negative integer n. 
We can then use the Lebesgue dominated convergence 
theorem to show that near z — >• +oo one has 



+oo 



(s)e~ 2s ~ z ds 



2n+l ' 



(B4) 



It now remains to specify the value of the constant Q 
for each component of D. The components x and y of 
D(co>) are equal and are denoted D xy . In particular, the 
terms connected to the two polarizations TE and TM are 
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denoted Dj£ E ' and D^ M ' respectively. They are given 
by 



D (TE) = °_ 

x v 4 

n(TM) = 3 
x v 4 



-2sz 



Im(/9i(s, J, e)) ds, 



—2sz 2 

e s 



Im(p2(^, S, e)) (is. 



(B5) 



For the z component D z of D(cj), we have 
Df E) = 0, 

4 TM) = t, / e" 2s " (1 + s 2 ) Im(p2(«, 5, c))d*. 
^ Jo 



(B6) 



The following asymptotic expansions near s = hold 

(B7) 



Im(pi(M,e)) ~-2sIm(/(£,e)) 
In 

where /(e, (5) 



Im(p 2 (s,^c)) ~ -2s Im(e /(£,€)), 

cot((5Ve ^ e then obtain for z —> oo 





D (TE) 


n(TM) 
J-^xy 


4™ ] 


5 = ^ ->. oo 


-3Im(/(e,<5)) 
85 2 


-9Im(e/(e,5)) 
165 4 


-3Im(e/(e,<5)) 
45 2 



Finally, we evaluate the behavior of the integrals with 
respect to S near and +oo. Denoting 5 = we have 
cot {8\/c - l) — >• ^t=y for (5 -» and cot (Sy/e - l) -» 
— z for (5 — >• oo, so that 



oo 



■ — »oo 



►0 



D 



xy 



85 



8z 2 S 



(TM) 



165 4 Re (v^l) 



9 



-9 



165 4 5 



Im 



(TM) 



Is Re 



4z- 



Az 2 5 



The integrals being continuous functions of the variables 
(z, 6) for z and 5 near -f oo, the limits z -f oo and 
5 — »> +oo commute. In conclusion, for z — » -f oo B(co>) 
is the only relevant contribution governing the interplay 
between the slab and environmental temperatures in the 
long-distance atomic dynamics. 

We now focus on the properties of D(cj) for small 
atom-slab distances z —> 0. As we will see, this evanes- 
cent contribution is divergent in this limit and is at the 



origin of the fact that the effective temperature at any 
frequency uo and thickness 5 tends to the one of the 
slab for small z. Indeed, the TE components are fi- 
nite because of the following behavior for s — » +oo: 
Im(pi(s, 5, e)) ~ Im(e)/4s 2 . The divergence comes in- 
stead from the TM components. In order to obtain the 
divergent TM small-distance behavior we use the follow- 
ing asymptotic theorem [42]: let h be a function such 
that h(s) has a polynomial behavior near s = -f oo, i.e. 
h(s) ~ H s n for some non-negative integer n. Then when 
z — >• 0, it holds that 



/ 

Jo 



e- 2zs h(s) ds~U 



2n+l gn+l ' 



(B8) 



The quantity Im(p2(s, S : z)) has expansions around s 
-f oo given by 



lm(p 2 (s,6,e)) ~ Im (^-py)- 
We then obtain 



(B9) 





xy/z 


n(TM) 


D (TM) 


z -> 


No divergence 


3Im(f=l) 


3Im(f=l) 


16z 3 


8« 3 



We observe that in this limit there is no dependence on 
slab thickness 8. 

The integrands being not continuous near (z,5) ~ 
(0, 0) the limits z — >• and (5 — >• do not commute (the 
notation 6 — » a, z — >• 6 means that first the limit z ^ b 
is taken, and then the limit 5 — » a is taken) and we have 





D (TM) 


n (TM) 
J-^xy 


J-^xy 


5 -> 0, z 


3 1 r 


3 1 r 




z 0, (5 


3 5 r 
64 5 3 i2 


3 <5 r 
128 5 3 i2 




S — » -f oo, z -> 


3 1 r 
8F Jl 


3 1 T 


!/n + °° Im(ri)ds 


z -» 0, J -> -f oo 


3 1 r 


3 1 r 

16 53 i l 


!Jn + °° Im(ri)ds 



where /i = Im (fqpi) an d ^2 = Im (j~T~^) ' 
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